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Abstract. We present a companion of the < 2 Myr young classical T Tauri star GQ Lup in the Lupus star forming 
region at 140 ±50 pc from imaging, astrometry, and spectroscopy. With direct K-band imaging using VLT/NAGO, 
we detected an object 6 mag fainter than GQ Lup located 0.7” west of it. Gompared to images obtained 2 to 
5 years earlier with Subaru/GIAO and HST/PC, this object shares the proper motion of GQ Lup by 5 and 7a, 
respectively, hence it is a co-moving companion. Its K — L' color is consistent with a spectral type early to mid L. 
Our NAGO K-band spectrum yields spectral type M9-L4 with H 2 O and GO absorption, consistent with the new 
GAIA-Dusty template spectrum for log g ~ 2 to 3 and Ten — 2000 K with ~ 2 Rjup radius at ~ 140 pc, hence 
few Jupiter masses. Using the theoretical models from Wuchterl & Tscharnuter (2003), Burrows et al. (1997), and 
Baraffe et al. (2002), the mass lies between 1 and 42 Jupiter masses. 
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1. Introduction: Direct detection of planets 

Direct detection of sub-stellar companions is difficult be¬ 
cause of the large dynamic range between the faint com¬ 
panion and the close-by, much brighter star. Due to con¬ 
traction, Myr young sub-stellar objects are brighter than 
Gyr older ones. Hence, young nearby stars would be the 
best targets for the direct detection of sub-stellar compan¬ 
ions. The mass of a detected sub-stellar companion can be 
estimated from the observed companion magnitude and 
the assumed or known age and distance of the primary 
star, using theoretical model calculations. Ghauvin et al. 
(2004) reported the detection of a faint object near the 
brown dwarf 2M1207, a potential member of TWA at ~ 8 
Myrs, whose infrared (IR) color and H-band spectrum are 
consistent with an L dwarf, hence ~ 5 Mjup (Baraffe et al. 
2002 models), if bound (astrometry missing). 

Here, we present evidence for a sub-stellar companion 
around GQ Lup located in the Lupus I cloud (Tachihara 
et al. 1996), a young K7eV-type classical T Tauri star with 
low extinction, but with mid- and far-IR excess (Hughes 
et al. 1994), i.e. a disk, and also both soft and hard X-ray 
emission (Krautter et al. 1997), rare for classical T Tauri 
stars (Neuhauser et al. 1995), so that we included this star 
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in our planet search programs by radial velocity and di¬ 
rect imaging since March 1999. As distance towards GQ 
Lup in the Lupus I cloud, we use 140 ± 50 pc (Wichmann 
et al. 1998, Neuhauser & Brandner 1998, Knude & Hog 
1998), the age of GQ Lup is 0.1 to 2 Myrs, depending 
on the set of models used (we work with 1 ± 1 Myr). 

2. Direct observations of a wide companion 

2.1. AO imaging detection and photometry 

We observed GQ Lup with Yepun (ESO-VLT-UT4) on 
Gerro Paranal, Chile, using the Adaptive Optics (AO) in¬ 
strument NAGO (Naos-Conica, Rousset et al. 2003) on 25 
June 2004 in visitor mode using the S13 camera (14” x 14" 
field) in the Kg-band. We took 200 co-adds x 0.347s x 27 
jitter positions. Data reduction was done with eclipse: 
dark subtraction, flat devision, shift±add. A companion 
candidate is found 0.7" west of GQ Lup (Fig. 1). By com¬ 
parison with the primary (K=7.096±0.020 mag, 2MASS), 
we got Kg = 13.1 ± 0.1 mag for the companion candidate 
with aperture photometry after subtracting the primary. 
K- and L’-band images of GQ Lup were obtained with 
the Subaru Coronagraphic Imager with Adaptive Optics 
(CIAO, Murakawa et al. 2004), in program o02312, re¬ 
trieved by us from the public Subaru archive SMOKA. 
We reduced the data in the above way. The candidate 
is detected in both K (as with NAGO) and L', where 
we obtain L' = 11.7 ± 0.3 mag by comparing the com- 
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Fig. 1. VLT-NACO Kg-band image of GQ Lup and its 6 
mag fainter companion candidate 0.7325 ± 0.0034" west. 


Table 1. Observing log, FWHM, and separations 


Tel./Instr. 

Obs. Date 

FWHM 

Separation 


HST/PG 

10 Apr 1999 

95 mas 

739 ±11 mas 

Subaru/CIAO 

17 Jul 2002 

150 mas 

736.5 ± 5.7 

mas 

VLT/NAGO 

25 Jun 2004 

68 mas 

732.5 ± 3.4 

mas 

VLT/NAGO 

25 Aug 2004 

74 mas 

731.4 ±4.2 

mas 

VLT/NAGO 

14 Sep 2004 

72 mas 

735.8 ± 3.7 

mas 


panion candidate with the primary.^ The primary has 
Ay = 0.4 ± 0.2 mag (Batalha et al. 2001), hence Ak = 
0.04 ± 0.02 and Ap' — 0.02 ± 0.01 mag (Rieke & Lebofsky 
1985). Applying this extinction to the companion, we de¬ 
rive {K — L')o ~ 1.4 ± 0.3 mag, consistent with L2-7 
(Stephens et al. 2001, Golimowski et al. 2004).^ At ~ 140 
pc, both the primary and the companion candidate are by 
~ 3 mag brighter than main-sequence K7V stars and Gyr 
old (main-sequence) early to mid L-dwarfs as far as the 
absolute magnitudes are concerned, consistent with both 
being young. 

GQ Lup and its companion candidate are also detected 
in archived images obtained with the Hubble Space 
Telescope (HST) Wide Field Planetary Gamera No. 2 
(WFPC2) in program SNAP 7387 (retrieved by us from 
the public MAST archive at STScI) in the filters F606W 
and F814W. Data reduction was performed for the HST 
data as in Neuhauser et al. (2002). Since the companion 
candidate is detected only marginally, we cannot obtain 
reliable optical colors. 


^ L= 6.05 ± 0.13 mag from Glass & Penston (1974) and 
Hughes et al. (1994); hence, L' ~ 6.05 mag; the difference 
between L and L' for K7 is marginal, within the 0.3 mag error 
given of the companion. 

^ The primary has a K-band excess due to its disk of ^ 2 
mag (from a K7V blackdody ht to its UBVRIJHK mags). 



Fig. 2. Observed separation between primary star GQ 
Lup A and companion candidate with HST/PC (left), 
Subaru/CIAO (middle) and VLT/NAGO (right, first the 
June 2004 image, then the Aug and Sept 2004 aquisition 
images, all consistent within Icr) compared to the expec¬ 
tation when bound (no change in separation as straight 
full line ±5 mas/yr for orbital motion as straight dotted 
lines) or when background (motion of GQ Lup A alone to 
the SW due to its proper motion ± its error with wobble 
due to expected parallactic motion of 7.1 mas). We can 
reject the background hypothesis by 4.9cr (NACO-CIAO) 
plus 7.2cr (NAGO-HST). 

2.2. Astrometry 

To study whether primary and companion candidate form 
a common proper motion pair, we compare separations 
and position angles PA at the different epochs. We use 
21.33 ± 0.02 mas/pixel as CIAO pixel scale (Fukagawa et 
al. 2003, 2004) and 45.545 ±0.005 mas/pixel for HST/PG 
(Holtzman et al. 1995). The positions of primary and com¬ 
panion are determined by Gaussian centering, after pri¬ 
mary PSF subtraction as far as the companion is con¬ 
cerned, and the errors given include the pixel scale error. 
The separations measured are given in Table 1 (Ict er¬ 
rors) including the aquisition images for the spectroscopy, 
obtained with the NACO S27 camera with 27.07 ±0.11 
mas/pixel. The astrometric precision is best in our own 
deep NAGO image (June 2004) obtained with the S13 
camera with 13.23 ± 0.05 mas/pixel. Scale and orienta¬ 
tion of S13 and S27 were obtained by us using a Ori AD 
(Hipparcos: 12.980" at 84.2°) and GJ 852 AB taking into 
account the orbital motion within them. The position an¬ 
gle of the detector was tilted by 0.14 ± 0.25° to the east. 
If the faint object would be a background object, it would 
have negligible proper motion and we would see a change 
in separation due to the proper and parallactic motion of 
GQ Lup A (Fig. 2). Its proper motion is /Tq, = —27 ± 3 
mas/yr and ns = —14 ± 3 mas/yr (Teixeira et al. 2000, 
Gamargo et al. 2003) at 140 ± 50 pc distance. The back¬ 
ground hypothesis can be rejected by 4.9 and 7.2a by com¬ 
paring NACO with CIAO and HST, respectively. The PA 
of the companion is 275.45 ± 0.30° in the NACO image in 





June 2004 (detector orientation calibrated); if the candi¬ 
date would be a non-moving background object, the PA 
should have been 269.89 ± 1.03° at the time of the HST 
epoch, but we observe 275.62 ±0.86° (detector orientation 
from fits file header, known to be stable and precise), i.e. 
4 . 3(7 deviant from the background hypothesis. Hence, the 
faint object is co-moving with GQ Lup A. The probability 
to find by chance a fore- or background M or L dwarf with 
the same proper motion near GQ Lup A is negligible. 

2.3. Spectroscopy of the companion 

To check whether the companion is indeed cool, we ob¬ 
tained a K-band spectrum with VLT/NAGO using the 
S54-4-3K mode from 1.79 to 2.57 ^m with resolution 700 
and nodding along the 172 mas slit. Two spectra were ob¬ 
tained on 25 Aug and 14 Sept 2004 with 20 and 40 spectra, 
respectively, exposed for 30s each and reduced in the nor¬ 
mal way: dark subtraction, flat devision, shift±add, and 
wavelength calibration (Fig. 3). 

We confirmed that the faint object was positioned 
well in the center of the slit in every individual spec¬ 
trum. Wavelength dependant refraction and Strehl-ratios 
through the narrow slit can change the slope of the spec¬ 
trum: The light-loss due to refraction, following Schubert 
& Walterscheid (2000) with the given airmasses, paral¬ 
lactic and positional angle, is 1.5 and 1.6% at the 1.9 
and 2.5/xm, respectively, compared to the flux at 2.35/im. 
Given the seeing and slit-width, the slit efficiency is 76.4% 
at 1.67/rm and 79.0% at 2.18/xm, so that the flux in the 
blue is reduced by 2.5%. In total, 2.0% of the flux is miss¬ 
ing in the blue compared to the middle of the spectrum 
and only 0.8% in the red. The slit-loss is partly compen¬ 
sated by higher Strehl-ratio, because the companion flux 
peaks in the middle of the K-band. 

The flux of the companion is calibrated once with the spec¬ 
trum of the K7eV primary star GQ Lup A observed simul¬ 
taneous and once with the telluric standard HD 159402 
(B3HI), observed in the same night. The large telluric 
absorption features below 2.06 pm and above 2.42 pm 
nevertheless cause a considerable amount of noise in this 
part of the spectrum. The Reid et al. (2001) K1 spectral 
index is 0.13 to 0.39 giving M9-L3. The McLean et al. 
(2003) H 2 O-D ratio is 0.67 to 0.89 indicating L2-7 (c.f. 
Cushing et al. 2005). The Nal doublet has Wa = 4.0 ±0.5 
A for the primary, but < 2 A for the companion (possi¬ 
bly partly telluric), which implies a spectral type ~ M9 or 
later (Cushing et al. 2005, Comeron et al. 2000) or ~ L2 or 
later (Gorlova et al. 2003); even for strongly reddened M 
stars, the Nal line should still be strong (Greene & Lada 
1996). The CO band head at 2.295/im is also present in 
M-dwarfs, increases slightly in strength for early L, but 
weakens again in later types. This band head is clearly 
visible with CO-index ~ 0.86, i.e. M 6 to LO (Gorlova et al. 
2003). The average of all the above estimates is L1.5±2.5, 
or M9 to L4. This is consistent with the dereddened K — L' 
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Fig. 3 . Our flux-calibrated K-band NACO spectrum of 
the companion of GQ Lup (bottom) compared to the 
young M 8 brown dwarfs in Cha I (top) from Comeron 
et al. (2000), L2 (2MASSW J0829066±145622, 2nd from 
top) from Reid et al. (2001), and a GAIA-Dusty template 
spectrum for 2000 K and log (/=2, which compares well 
with the companion. 

color; hence no evidence for additional extinction due to, 
e.g., a disk around the companion. 

3. Mass determination and discussion 

Our companion, bound to a 1 ± 1 Myr young star, may 
be the youngest and lowest-mass companion ever imaged. 
Hence, it is difficult to compare it to field (old) L-dwarfs. 
From the K-band CO-index being 0.862 ± 0.035 we ob¬ 
tain the gravity log g = 2.52 ±0.77, a slight extrapolation 
from late-M to early-L (Gorlova et al. 2003), hence uncer¬ 
tain. Since neither spectral type nor gravity are well con¬ 
strained, the temperature Teff is only weakly constrained 
to ~ 1600 to 2500 K (or 2050 ±450 K) from Gorlova et al. 
(2003) and references therein as well as Reid et al. (1999), 
Stephens et al. (2001), Burgasser et al. (2002), Nakajima 
et al. (2004), and Golimowski et al. (2004). 

We compared our observed spectrum with the theoretical 
template spectra from the GAIA-Dusty model (Brott & 
Hauschildt, in prep.), updated from Allard et al. 2001, 
with improved molecular dissociation constants, more 
dust species with opacities, spherical symmetry, and a 
mixing length parameter 1.5-Hp; we tried Tea = 2000 and 
2900 K and log g = 0, 2, and 4. A good fit is obtained 
only for 2000 K and log g = 2 (Fig. 3). Because the fit 
for log g = 0 is much worse than for log 5 = 4, where 
the continuum at 2.22 to 2.3 pm and the depth of the CO 
lines are not reproduced, and because of the log g from 
the CO-index (see above), we conclude log (/ ^ 2 to 3. 
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Fig. 4. H-R diagram with tracks from Wuchterl & 
Tscharnuter (2003) plus more tracks calculated by us 
(from top to bottom) for 1 and 0.7 Mq as well as 13, 
5, 2, 1, and 0.5 Mjup. The isochrones (dashed) for 1.10 
and 1.65 Myr connect the end of the 0.7 M 0 track (star 
A) with the planetary tracks. Components A and B (Itr 
errors) are co-eval at ~ 1.1 Ma. 


The observed flux can be reproduced for an object with 
~ 2 Rjup radius at 140 pc. 

With B.C.k = 3.3±0.1 mag for M9-L4 and K-L’~ 1.4±0.3 
(Golimowski et al. 2004), Kg = 13.1 ± 0.2 mag for the 
companion, and 140 ± 50 pc distance, the luminosity is 
\og{L/L q) = —2.37±0.41 for the companion. We can now 
plot it into an H-R diagram to compare with models. From 
Burrows et al. (1997) Fig. 7, the companion has ~ 12 to 32 
Mjup (from L and age), but only ~ 3 to 9 Mjup from Figs. 9 
& 10 (from T and age). It is similar in Baraffe et al. (2002) 
Fig. 2: The companion has ~ 3 to 16 Mjup from T and age, 
but ~ 11 to 42 Mjup from L and age; ^ 12 to 42 Mjup from 
Mk and age (see also perso.ens-lyon.fr/isabelle.baraffe). 
Burrows et al. and Baraffe et al. start with an assumed 
internal structure without collapse, so that their models 
are uncertain up to a few Myrs (Baraffe et al. 2002). 
Wuchterl & Tscharnuter (2003) include the initial col¬ 
lapse (Fig. 4): Tracks for masses of 1 to 0.013 Mq are 
radiation fluid-dynamical calculations of the collapse of 
initially marginally unstable Bonnor-Ebert-spheres; plan¬ 
etary tracks are models in the framework of the nucleated 
instability hypothesis (Wuchterl et al. 2000; Wuchterl, 
in preparation). All ages are counted from first photo¬ 
sphere formation. The age-offset visible on the 1.10 Myr 
isochrone, at the transition from the 0.013 Mq brown 
dwarf collapse to the 5 Mjup planetary track is due to the 
time-offset caused by the planetary core formation (< 10® 
yrs). The companion has a mass of ~ 1 to 2 Mjup (lu) and 
is co-eval with the star (~ 1.1 Myr). From this model (Fig. 
4), we expect log g ~ 2.4 and ~ 1.8 Rjup radius, consis¬ 


tent with our spectrum and the best fitting GAIA-Dusty 
model spectrum (Fig. 3). 

The most critical point in the mass determination of the 
companion (candidates) of GQ Lup and 2M1207 are the 
models, which may be off by an unknown factor for low 
ages (few Myrs); they need to be calibrated, before the 
mass of such companions can be determined confidently. 
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